PROCEEDINGS 

°F SCIENCE 



Non-perturbative renormalisation of four-fermion 
operators in A^f = 2 QCD* 



o 
o 

(N 

> 

o 
O 

(N 



(N 
> 
(N 

oo 



O 



Hlpha 

Collaboration 

Petros Dimopoulos, Gregorio HerdoizaJ Anastassios Vladikas 

INFN, Sezione di Roma "Tor Vergata" and Dipartimento di Fisica, Universitd di Roma "Tor 
Vergata " 

Via delta Ricerca Scientifica 1, 1-00133 Rome, Italy 
email: petros . dimopoulos@roma2 . inf n . it , 

gregorio . herdoiza@roma2 .infn.it, tassos. vladikas@roma2 . inf n . it 

Filippo Palombi, Mauro Papinutto 

CERN, Physics Department, Theory Division 
CH-1211 Geneva 23, Switzerland 
email: 



filippo . palombigcern . ch, mauro . papinutto@cern . ch 



Carlos Pena* 

Departamento de Fisica Teorica and Instituto de Fisica Tedrica UAM/CSIC 

Facultad de Ciencias, Universidad Autdnoma de Madrid, Cantoblanco E-28049 Madrid, Spain 



email: carlos.pena@uam.es 



Hartmut Wittig 

InstitutfUr Kernphysik, University of Mainz 
D- 5 5099 Mainz, Germany 



email: »fittig@kph . uni-mainz . de 



We present results for the non-perturbative renormalisation of four-fermion operators with two 
flavours of dynamical quarks. We consider both fully relativistic left current-left current opera- 
tors, and a full basis for AZ? = 2 operators with static heavy quarks. The renormalisation group 
running of the operators to high energy scales is computed in the continuum limit for a family 
of Schrodinger Functional renormalisation schemes, via standard finite size scaling techniques. 
The total renormalisation factors relating renormalisation group invariant to bare operators are 
computed for a choice of lattice regularisations. 

The XXV International Symposium on Lattice Field Theory 
July 30-4 August 2007 
Regensburg, Germany 



*CERN-PH-TH/2007-178, FTUAM-07-14, IFT-UAM-CSIC-07-45. 

'Address after Oct 1st: DESY-Zeuthen, Platanenallee 6, D-15738 Zeuthen, Germany. 

''"Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



NPR offour-fermion operators in Nf = 2 QCD 



Carlos Pena 



1. Introduction 

Hadronic matrix elements (HMEs) of four-fermion operators have long been essential input 
quantities for Flavour Physics. Reliable estimates of a number of HMEs are crucial in the study 
of CP violation via CKM unitarity triangle analyses, or of such striking experimental findings as 
the enhancement of hadronic decay amplitudes by long-distance effects (as e.g. in the A/ = 1/2 
rule). The only known technique to compute HMEs from first principles is lattice QCD. However, 
lattice QCD results have long been hampered by the difficulty to eliminate a number of systematic 
uncertainties. Most notably, the high cost of including dynamical quark effects in lattice QCD sim- 
ulations has enforced for many years either the quenched approximation, or the use of dynamical 
quark masses far too heavy to allow for a well-controlled extrapolation to the physical regime. In 
some cases, e.g. the computation of the kaon bag parameter Bk, quenching effects are indeed the 
last remaining uncontrolled systematic uncertainty [|T]]. 

As techniques for the simulation of light dynamical quarks have witnessed dramatic progress 
in the last few years (see e.g. [Q]), it becomes increasingly important to bring to this environment 
the techniques to control other sources of uncertainty, in order to aim at precision computations 
of physical quantities. In the context of HMEs, one of the most prominent examples is non- 
perturbative renormalisation (NPR) (see e.g. [||]). The use of finite-size scaling techniques has 
allowed to control fully both the renormalisation group (RG) running and the matching of lattice 
to renormalised observables in the quenched approximation for a broad class of four-fermion oper- 
ators [Q, ^, ^]. The aim of the present work is to extend these results to N{ = 2 QCD. In particular, 
we will discuss 1 . the RG running of left current-left current relativistic four-fermion operators, 2. 
the RG running of all AB = 2 operators with static heavy quarks, and 3. the matching of the above 
operators to renormalised continuum operators for some particular choices of the regularisation. 
Immediate applications, as we will point out later, arise in the computation of the bag parameters 
Bk and Bb- Preliminary results had been presented at last year's conference [^. 

2. Definitions and setup 

2.1 Renormalisation of four-fermion operators 

We will consider two different classes of operators: 

Or,.r, W = \ [{¥i{^)riV2{x)) mx)r2Y4{x)) ± (2 ^ 4)] , (2.1) 
^r,,r, W = \ [{%{x)riY2{x)) (%(x)r2i/A4(x)) ± (2 ^ 4)] . (2.2) 

In the above expressions xjfk is a relativistic quark field with flavour index k, Xffh h static (anti)quark 
fields, F/ are spin matrices, and the parentheses indicate spin-colour traces. All the fields are 
interpreted to be in the valence sector of the theory. This formalism of distinct quark flavours 
will allow us to isolate scale-dependent logarithmic divergences from eventual mixing with lower- 
dimensional operators that may appear for specific choices of quark masses and/or flavour content. 

The above operators mix under renormalisation as determined by the symmetries of the regu- 
larised theory. If we restrict ourselves to the parity-odd sector, complete bases of operators in the 
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relativistic and static cases are given by 

(1^ (= -fo^ (1^ (1^ (1^ (1^ \ ■ (= I ff^ ff^ ff^ iff^ \ (? '^^ 

^ \'-^VA+AV) '-^VA-AV) "-^SP-PS) "-^SP+PS) '-^TT' / ' "^k ^ \ ^ VA+ AV ) ^ VA- AV ) SP-PS > ^ SP+PS 5 / ' V^--"/ 

respectively, in standard self-explanatory notation for the choice of spin matrices F/. A full analysis 
of the renormalisation of these operator bases with relativistic Wilson fermions has been performed 
in [§, One particular conclusion is that, contrary to the parity-even case, discrete symmetries 
protect all the above operators from extra mixings under renormalisation due to the breaking of 
chiral symmetry. Recall that the RG of these operators and of their parity-even partners is identi- 
cal, as in the continuum limit (CL) chiral symmetry holds. On the other hand, the connection to 
observables involving matrix elements of parity-even operators is non-trivial. 
From now on, we will consider the subset of operators 

Ql : '^k ^ ~1~ 4c^2 5=^3 ? 2c^4 . (2.4) 

All these operators renormalise multiplicatively — i.e., given an operator O G {2f ,=S^^} the cor- 
responding operator insertion in any on-shell renormalised correlation function is given by 

0^{x;ijl) = limZ(^o,aAi)C(x;^o), (2.5) 

where gQ,a we, the bare lattice coupling and the lattice spacing, respectively. The RG running of 
the operator is controlled by the anomalous dimension 7, defined by the Callan-Symanzik equation 

^j^Or{x-^) = 7(g(Ai)) ORix;n) , (2.6) 

which is supplemented by the corresponding Callan-Symanzik equation for the renormalised cou- 
pling 

At^f(M)=iS(f(At)). (2.7) 

In mass-independent renormalisation schemes, the beta function and all anomalous dimensions do 
indeed depend only on the renormalised coupling g. They admit perturbative expansions of the 
form 

j8 (g) -g' {bo + b,g^ + b2g^ + ...); rig) -g^ (70 + Tig^ + r2g' + ■■■), (2.8) 
in which the coefficients bo,bi,Yo are renormalisation scheme-independent. Upon formal integra- 



tion of Eq. (2.6), one is left with the renormalisation group invariant (RGI) operator insertion 



/ /■«"". /r(s) »\1 „„, 



471 

while the RG evolution between two scales /ii , /X2 is given by the operator 



t/(At2,Mi)=exp dg^ =lim-p^. (2.10) 
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2.2 Schrodinger Functional renormalisation schemes 



Eq. ( |2. 10| ) is the starting point to compute non-perturbatively the RG evolution of composite 
operators. To that purpose we introduce a family of Schrodinger Functional (SF) renormalisation 
schemes. The latter are defined by regularising the theory on a symmetric lattice of physical size 
with SF boundary conditions (see e.g. [|T^ for an introduction to the SF setup). The renormalisation 
scale is set to be the infrared cutoff, i.e. jj. = l/L. Renormalisation conditions for relativistic 
operators have the form 

F{xo) F{xo) 



?o,alJ.) 



& & 



(2.11) 

tree level 



and are imposed in the chiral limit. In the above expression, F is a four-point correlation function 
of the form 

F(xo) = ■^(^2l[rA]^45[rB](2,^(x)^^3[rc]) , (2.12) 

where ^[T] are bilinear interpolating fields living on the time boundaries, and is a suitable 
boundary-to-boundary correlation function that divides out the ultraviolet divergences associated to 



these bilinears. Similar renormalisation conditions to Eq. (2. 1 1) are set up for static-light operators, 
with flavours 1 and 3 substituted by h and h. Full details are provided in |Q, ^, ^. For now it is just 
important to mention that the renormalisation scheme is fully determined by fixing the parameters 



involved in the SF boundary conditions; the point xq at which Eq. ( [2. 1 1[ ) is imposed; the Dirac 
matrices Ta^b.c entering boundary bilinears ^; and the normalisation factor 0. Specific schemes 
have been introduced in [Q, ||, Here we will concentrate in the cases which have been found 
to be best behaved in the quenched study, namely scheme 1 for and scheme 8 for Qy in the 
notation of [Q, ^, and the reference schemes for static-light operators defined in |^. 

A crucial observation is that all the above renormalisation schemes are mass-independent by 
construction, and the resulting renormalisation factors are flavour-blind. It then follows that they 
can be used to remove the logarithmic divergences from any four-fermion operator with the consid- 
ered structure, irrespective of its specific flavour content, once eventual subtractions due to mixing 
with equal- or lower-dimension operators have been properly performed. 

2.3 Step-scaling functions 

The basic objects to study the RG evolution of composite operators non-perturbatively are the 
step-scaling functions (SSFs) 



Z{go,a/L) 



(2.13) 

gHUL)=u 



which can be computed at several values of the lattice spacing for fixed physical size (inverse 
renormalisation scale) L. The corresponding values of j3 are indeed fixed by requiring that the 
renormalised SF coupling g^, and hence L, are kept constant. It is then possible CL extrapolation 

a{u)^l\mL{u,a/L) = U {{2L)-\L-')\_,^^^^^^^ . (2.14) 



At vanishing external momenta, thiere are 5 possible nontrivial choices that preserve cubic symmetry. 
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Once a{u) is known for several different values of the squared gauge coupling u, it is possible 
to reconstruct the RG evolution factor U (/ihadiMpt) between two extreme scales /Xhad> in the range 
of a few hundred MeV, and /Xpt in the high-energy regime. This in turn allows to compute the RGI 
operator of Eq. ( ^!^ ) in a way free from large uncontrolled systematic uncertainties. It is enough to 
consider the exponential on the rhs of Eq. (|2.9|) evaluated at n = /ihad> and split it as 



The second factor on the rhs is known non-perturbatively, while the first factor can be safely com- 
puted at NLO in perturbation theory, provided the scale |J,p^ is high enough so as to render NNLO 
effects negligible. 

3. Non-perturbative computation of the RG running 

SSFs have been computed using the non-perturbatively 0(a) improved Wilson action, and 
a HYP2 action for static quarks, at six different values of the SF coupling, corresponding to six 
different physical lattice lengths L. For each volume we have simulated at three different values 
of the lattice spacing, corresponding to lattices with L/a = 6, S, 12 (respectively L/a = \2, 16,24) 
for the computation of Z(L) (resp. Z(2L)). We used the A'f = 2 configurations generated by the 
ALPHA Collaboration for the determination of the RG running of the quark mass [|n|]. All the 
technical details concerning the dynamical simulations are discussed in the mentioned work. 

As we do not implement full 0(a) improvement for four-fermion operators, the only linear 
cutoff effects that are removed from L{go,a/L) are those cancelled by the SW term in the fermion 
action. Therefore, we expect SSFs to approach the CL linearly in a/L. In practice, it is often 
observed that the data corresponding to L/a = 8, 12 are compatible within errors, whereas the 
L/a = 6 datum, that is expected to bear the largest cutoff effect, is off. This suggests that a weighted 



average of the results for the two finest lattices, as considered in []llp, may yield a good estimate 
of the CL value. However, the lack of at least one extra value of a/L closer to the continuum, 
that would allow a more precise control of the systematics, has led us to conservatively adopt 
linear CL extrapolations involving all the data. It is worth remarking, though, that linear fits and 
weighted averages lead to compatible results within one standard deviation in most cases, as can 
be seen in Fig. |l|. The latter illustrates the extrapolations at all values of the coupling for two 
selected operators. Finally, let us mention that autocorrelation times, which are included in the 
error estimate, increase towards the CL, leading to amplified errors in the finest lattices. 

The resulting SSFs a{u) have been fitted to a polynomial form. For definiteness, we will 
provide results for a fit to a{u) = 1 + siu + S2U^ + st,u^ , where si is fixed at the value predicted 
by LO perturbation theory and S2,S3 are left as free parameters. Once this continuous form of 
the SSF has been obtained, it is possible to compute the relation between the RGI operators and 
the renormalised operators at the low-energy scale /Xhad = ^max' defined by g^(Lmax) = 4.61, as 
explained e.g. in [^, |^. This scale is chosen such that the renormalisation constant Z{go,aiJ.had) 
can be computed on accessible lattices in ranges of values of go commonly used in large volume 
simulations. The results for the operators under investigation are reported in Table |l[ Note that 
typical relative errors reach the 5% ballpark, which may result in a sizeable error in HMEs coming 
from renormahsation alone. 
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Figure 1: CL extrapolation of the SSFs for (left) and (right) for one particular choice of the 
renormaUsation scheme (all boundary matrices set to 75, spatial boundary conditions set by = 0.5, a = 
in the renormalisation condition for (see ^ ^] for details). The renormalised coupling increases from 
top to bottom and from left to right. Blue discontinued lines and the blue point at a/L — Q correspond to 
weighted averages of the L/a = 8, 12 data, red lines and the red a/L = cross to linear extrapolations in a/L 
of the three data. 





1.0 2.0 3.0 



Figure 2: SSFs for (left) and (right) in the CL in the same schemes as considered in Figure 1. Full 
red points are Nf = 2 results, open blue points are quenched results. The red dashed (blue dotted) line is the 
NLO A^f = 2 (A^f = 0) perturbative result. 

4. Connection to hadronic observables 



RGI operator insertions can be related to bare operator insertions via a total renormalisation 
factor Z{go), defined as 



'dgfM_Jl 



(4.1) 



This is enough to remove all ultraviolet divergences, once eventual renormalisation scale-independent 
mixing with operators of dimension d <6 peculiar to the specific flavour structure under consid- 
eration has been taken into account via suitable subtractions. The details of the mixing depend on 
the regularisation in which bare correlation functions are computed, as does the relation between 



the latter and physical observables. For instance, in [ ]12| , |13[ ] it has been explained how to extract 
the bag parameters Bk and Bg (the latter in the static Umit for the b quark) directly from three-point 
functions involving the operators and ^^2^ using Wilson actions with suitable twisted mass 
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operator 



ratio 



operator 



ratio 



1.201(66) 
0.554(21) 




0.724(34) 
0.647(32) 
0.539(18) 
0.796(20) 



Table 1: Ratios of RGI to renormalised at jihad operator insertions for the various operators in the reference 
renormaUsation schemes mentioned in the text. 

terms. The computation of the RGI renormaUsation factors Z{go) at a number of values of the bare 
coupling with an 0(a) improved Wilson action is under way and close to completion. 

5. Conclusions 

We have presented a fully non-perturbative computation of the RG running of a wide class 
of four-fermion operators in A^f = 2 QCD. These results, together with the matching to specific 
hadronic schemes, is a basic building block of any Nf = 2 computation of such quantities as Bk 
and Bb that aims at eliminating systematic uncertainties related to renormaUsation. On the other 
hand, the precision of the results sets a potentially unsatisfactory lower bound for the final error 
on weak matrix elements. Future refinement, e.g. by adding a finest lattice to our continuum limit 
extrapolations, can be hence desirable. These issues will be discussed in detail in our forthcoming 
pubUcation of the definitive results. 
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